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ABSTRAGT 

Post-starburst (or “E-l-A”) galaxies are characterized by low Ha emission and strong Balmer absorp¬ 
tion, suggesting a recent starburst, but little current star formation. Although many of these galaxies 
show evidence of recent mergers, the mechanism for ending the starburst is not yet understood. To 
study the fate of the molecular gas, we search for GO (1-0) and (2-1) emission with the IRAM 30m 
and SMT 10m telescopes in 32 nearby (0.01 < z < 0.12) post-starburst galaxies drawn from the Sloan 
Digital Sky Survey. We detect GO in 17 (53%). Using GO as a tracer for molecular hydrogen, and a 
Galactic conversion factor, we obtain molecular gas masses of M{H 2 ) = 1O®'®-1O® ®M0 and molecular 
gas mass to stellar mass fractions of ^ 10“^-10“°'^, comparable to those of star-forming galaxies. 

The large amounts of molecular gas rule out complete gas consumption, expulsion, or starvation as 
the primary mechanism that ends the starburst in these galaxies. The upper limits on M{H 2 ) for the 
15 undetected galaxies range from IO^'^Mq to IO^'^Mq, with the median more consistent with early- 
type galaxies than with star-forming galaxies. Upper limits on the post-starburst star formation rates 
(SFRs) are lower by ^ lOx than for star-forming galaxies with the same M{H 2 ). We also compare 
the molecular gas surface densities (Shs) to upper limits on the SFR surface densities (Ssfr), finding 
a significant offset, with lower User for a given Eh 2 than is typical for star-forming galaxies. This 
offset from the Kennicutt-Schmidt relation suggests that post-starbursts have lower star formation 
efficiency, a low GO-to-H2 conversion factor characteristic of ULlRGs, and/or a bottom-heavy initial 
mass function, although uncertainties in the rate and distribution of current star formation remain. 
Subject headings: galaxies: formation — radio lines: galaxies 


1. INTRODUCTION 

Post-starburst (or “E-l-A”) galaxies show signs of being 
caught in the middle of a dramatic, but brief, stage in 
their evolution. Emission line indicators suggest little-to- 
no current star formation, but strong Balmer absorption 
lines indicate a population of A stars that formed in a 
substantial b urst of star formation before a sudden stop 
^ I G yr ago ((Dressier fc Gunnlll983t iCouch fc Sharplea 

[MS). 

Post-starburst galaxies are likely in transition between 
star-forming gas-rich disk galaxies and passively evolv¬ 
ing gas-poor early types. Their disturbed morpholo¬ 
gies indicate that many are post- merger, and mos t have 
spheroid-dominated kinematics (IZabludoff et al.l II996I : 
[Norton et al.ll200ll : lYang et all 12004 1200811 . Many have 
blue cores, which can fade int o the color gradients ob¬ 
served in early t y pe ga laxies ([Yamauchi fc Gotol 120051 : 
lYang et all 120061 1200^ . and many lie in the “Kreen 
valley ” of the color magnitude diagram (jWong et al.l 
l2Q12li . Although only ^ 0.2% of local galaxies are post- 
starbursts, the short duration of this phase suggests 
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that ^ 40% of galaxies could have pas sed through it 
((Zabludoff et al.lfl996t ISnvder et al.ll20li[l . 

A critical part of galaxy evolution is the end, or pos¬ 
sible “quenching,” of star formation. As transitional ob¬ 
jects, post-starburst galaxies serve as a unique laboratory 
for understanding the processes that drive this cessation. 
Explanations for the end of the starburst fall into two 
general categories: elimination of the molecular gas or 
suppression of star formation. 

One possibility is that the starbu rst uses up the 
dense molecular cloud s in forming stars (iKennicuttK 1 9981 : 
iGao fc SolomonI [20041 . Molecular gas could also be re¬ 
move d from the galaxy in outflows ((Narayanan et ahl 
I2008I) . Evidence of LINER activity and la r ge outflows 
are s e en in post-starbursts ((Yan et al.1120061 : (Yang et aP 
l2006t [iLemonti et al.l 120071) . a nd AGN are observe d to 
eject molecular gas in outflows ((Feruglio et al.|[20T^ . al¬ 
though the driver of the outflows in post-starb ursts may 
be du e to star formation, not AGN activity ((Sell et al.1 
120141 . Some environmental effects, such as starv ation 
(e.g., iLarson et al.1119801: iBoselli fc: Gavazzil 1200611 . are 
thought to eliminate molecular gas reservoirs in galaxies. 
The molecular gas mass is several orders of magnitude 
lower in early types than in late types ((e.g.. lYonng et al.1 
120111: ICrocker et ffi]l2011h . If post-starbursts are becom¬ 
ing early types, they must lose or repurpose most of their 
gas. 

Feedback mechanisms could be responsible for sup¬ 
pressing star formation, resulting in the end of the star- 
burst. Molecular gas heating and suppressed star for- 
matio n efficiency have bee n claimed in galaxies with 
AGN ((Nesvadba et al.l[2010ll . resulting in higher observed 




































































2 


molecular gas sur face densities th an the Keirnicutt- 
Schmidt relation (|Kennicutt] I1998D would predict for 
their star formation rate (SFR) densities. Observations 
of cold gas in early type galaxies with AGN and re¬ 
cent bursts of star formation reveal little molecular gas 
(< 10® Mq), which declines st eeply with the age since the 
last period of star formation (iSchawinski et al1l2009f) . In 
our sample, the timescales necessary for outflows to expel 
the molecular gas from the galaxy are less than the time 
elapsed since the starburst ended (about 0.3-1 Gyr), so 
if AGN feedback has signihcantly reduced the molecu¬ 
lar gas reservoirs, we should observe the galaxies in their 
depleted state. 

A lower star formation efficiency is sugge sted in gas- 
rich, fast-rotating early type galaxies by iDavis et al.l 
(|2014f) . who observe lower SFR surface densities than 
the molecular gas surface densities would predict by a 
factor of ~ 2.5. These authors favor dynamical methods 
of lowering star formation efficiency in t his sample of 
galax ies, such as morphological quenching (|Martig et al.l 
1200911 . where the gravitational stability of the gas pre¬ 
vents it from collapsing and forming stars. Although 
post-starburst galaxies are likely to evolve into early 
ty pes, it is not clear that the gas-rich sample studied 
by I Davis et al.l (12014) are on the same evolutionary se¬ 
quence as post-starbursts. 

We aim to test these explanations for the starbursts’ 
end by constraining the properties of molecular gas 
within post-starbursts. Reservoirs of HI have been 
observed in post - starburst galaxies (|Chang et al.l I200ll : 
iBuvle et al.ll2006l : iZwaan et al.ll2013ll . In six of the eleven 
post-starbursts targeted in these samples, HI 21 cm emis¬ 
sion is detected, with atomic gas to stellar mass fractions 
typically between those of early and late type galaxies. 
However, HI is not a good tracer of star formation fuel 
(|Kennicutt. Jr. et al.ll2007n . and we must look at molec¬ 
ular gas signatures to understand the starbursts’ end. 

Detailed GO maps have been measure d for only a 
hand f ul of local post-sta rburst galaxies (|Kohno et al.l 
120021 : lAlatalo et al.l [2M^ . Even then, the two galax¬ 
ies studied, NGG5195 and NGG1266, are not universally 
agreed-upon as post-starbursts due to their Ha emis¬ 
sion. The molecular gas in these galaxies is centrally con¬ 
centrated, reaching starburst-like gas surface densities. 
Their kinematics l ed these authors t o suggest morpho¬ 
logical quenching (|Martig et al.l [200^ . where the gravi¬ 
tational stability of the gas prevents it from collapsing 
and forming stars. There is a need for a survey of the 
molecular gas content in a representative sample of post- 
starburst galaxies. 

We set out here to determine how much molecular gas 
remains in a sample of 32 post-starburst gal axies drawn 
from the Sloan Digital Sky Survey (SDSS, lYork et ^ 
[ 2003 ), and to determine whether the molecular gas den¬ 
sities are consistent with the small or negligible levels of 
current star formation. We observe the CO (1-0) and CO 
(2-1) lines with the IRAM 30m telescope, and observe 
a subset of 13 galaxies in CO (2-1) with the SMT 10m 
telescope. By assuming that the CO traces H 2 , we test 
whether the cessation of star formation was due to a lack 
of molecular gas, or to the gas being consumed by the 
burst, expelled in outflows, o r prevented from en tering 
the galaxy (starvation of HI (|Larson et al.l flMOri '). By 
comparing to the molecular gas vs SFR surface density 


relation for other galaxies, we will be able to determine if 
the star formation efficiency in post-starbursts is reduced 
by either gas heating, morphological quenching, or some 
other mechanism. 

We discuss our sample and observations in §2. Mea¬ 
surements of molecular gas masses and comparisons to 
the SFRs are presented in §3. We test these results and 
consider their implications for galaxy evolution in §4, pre¬ 
senting our conclusions in §5. When needed, we assume 
a cosmology of = 0.3, Ha = 0.7, and h = 0.7. 

2. OBSERVATIONS AND DATA ANALYSIS 
2.1. Sample Selection 

Our parent sample is dr awn from the SDS S main 
galaxy spectroscopic sample ([Strauss et al.ll200'^ , which 
is selected to have a limiting magnitude of r < 17.77 
mag. The initial sampl e was selected from the SDSS DR7 
([Abazaiian et al.ll2OCi90 . using the line fluxes and i ndices 
from the MPA-JHU catalogs (| Aihara et al.l [20Tl[l . We 
exclude galaxies with z < 0.01 to eliminate those that 
are very large on the sky relative to the 3" diameter of 
the SDSS fibers. We also exclude galaxies with unreliable 
Q Ha equivalent widths (EW), or median signal-to-noise 
(S/N) values of less than 10 per pixel. These cuts en¬ 
sure that the line index measurements are reliable. Our 
final parent sample from DRIO is composed of 595,268 
galaxies. 

We select post-starburst galaxies from our parent sam¬ 
ple by identifying galaxies with strong stellar Balmer ab¬ 
sorption lines signifying a recent (;< Gyr) starburst but 
little nebular emission indicative of on-going star forma¬ 
tion. We use the Lick HJ index to characterize the stellar 
Balmer absorption. We require HJa — o'(H(5a) > 4A, 
where (t(H(5a) is the measurement error of the Hi5a in¬ 
dex. We ensure that the galaxies have little on-going star 
formation by requiring Ha EW < 3 A in the rest frame. 
These selection criteria result in a sub-sample of 1207 
galaxies from the parent sample (0.20%). 

We have chosen two sub-samples for HST, Spitzer, and 
Herschel imaging, which we follow-up here. 15 galax¬ 
ies designated “S” throughout were selected to repre¬ 
sent a variety of ages since the end of the burst and 
based on their projected 8^m flux from SDSS spectra 
and serendipitous Spitzer observations. Galaxies with 
nearby companions and large [OHI] equivalent widths in¬ 
dicative of AGN activity were excluded. The post-burst 
ages are determined by fitting stellar population synthe¬ 
sis (SPS) models to the galaxy spectrum, assuming a 
combination of a young and old single burst stellar pop¬ 
ulations (French et al. in prep). 17 galaxies designated 
“H” throughout were selected from their bright WISE 
12pm fluxes and again for a range of post-burst ages (al¬ 
though without the [OHI] equivalent width cut). More 
details on the “H” and “S” selection processes are avail¬ 
able in Smercina et al. (in prep). The effect of these 
selection criteria on properties of the resultant sample 
is studied in EH Basic parameters of this sample are 
listed in Table 1. 

2.2. IRAM 30m CO Observations 

Observations were carried out with the IRAM 30m 
telescope over two observing campaigns in January 2012 

^ We require h_alpha_eqw_err > -1 
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(project ID: 218-11) and in August - September 2012 
(ID: 074-12). We use the Eight Mixer Receiver (EMIR) 
to observe both CO(l-O) and CO(2-l) lines (rest fre¬ 
quency: 115.271 and 230.538 GHz). For each target, we 
tuned the 3mm band (E090) and 1.3mm band (E230) re¬ 
ceivers to the redshifted GO(l-O) and 00(2-1) frequen¬ 
cies, L^ohs = 103.5 - 113.5 GHz and 207.1 - 227.1 GHz, re¬ 
spectively. EMIR provides a bandwidth of 4 GHz in dual 
polarization corresponding to ^11000 and 5500 kms“^ 
for CO(l-O) and CO(2-l) lines, respectively. The Wide¬ 
band Line Multiple Autocorrelator (WILMA) was used 
as the backend, with a resolution of 2 MHz corresponding 
to ^5 kms“^ in the 3mm band. Data were taken with a 
wobbler-switching mode with a frequency of 0.5 Hz or 
IHz with a throw distance of 120" in azimuth. The 
weather varied significantly: the precipitable water va¬ 
por (PWV) ranged from 1mm to 10mm with medians 
of 3mm (winter) and 6.8mm (summer). Galibration was 
performed every 15 min with standard hot/cold load ab¬ 
sorbers. The pointing was checked every 2 hours and was 
found to be stable within 3". The FWHMs of beam are 
Ri22" and 11" for CO(l-O) and GO(2-l) lines, respec¬ 
tively. 

We reduced the data with CLASS within the GILDAS 
software packag(0 and IDL routines. We use the velocity 
intervals [-1200, -400] and [400, 1200] kms“^ to fit first 
order polynomials for baseline subtraction. The spectra 
are coadded weighted by the rms noise of each scan. The 
on-source time (Tqn) ranges from 12 to 100 min depend¬ 
ing on the strength of the line toward the targets. If 
the source was not detected within 3 hours at the tele¬ 
scope (Ton ~ Ihr), we moved on to the next target. The 
resulting rms noise per 5kms“^ bin are 1.1 - 4.4 mK 
and 1.9 - 9.7mK for the GO(1“0) and GO(2-l) obser¬ 
vations, respectively (Ta* scale). The conversion factors 
from K {Ta scale) to Jy at our observed frequencies are 
'^7.7JyK“^ and ~6.0JyK“^ for the 1.3mm and 3mm 
bands, respectively. We summarize the IRAM 30m GO 
observations in the Tables [2] and [3] 

2 . 3 . SMT CO Observations 

Observations at the SMT 10-m telescope were per¬ 
formed over four runs in May 2011, February 2012, De¬ 
cember 2012, and February 2013. We used the 1mm 
ALMA Band 6 dual polarization sideband separating SIS 
receiver and IMHz filterbank to measure the GO(2-l) 
230.5 GHz (redshifted to 207.1 - 227.1 GHz for our sam¬ 
ple) line for 13 post-starburst galaxies. The beam size 
of the SMT for this line is ~33". Beam switching was 
done with the secondary at 2.5 Hz and a throw of 120 ". 
Galibration using a hot load and the standard chopper 
wheel method was performed every 6 minutes. Calibra¬ 
tion using a cold load was performed at every tuning. 

To reduce the data, we again use CLASS. The main 
beam efficiency rimb is calculated using Jupiter in 
each polarization. We subtract a first-order polyno¬ 
mial baseline from the spectrum using data between 
[—500, 500] kms“^, excluding the central region of [—300, 
300] kms“^. The spectra are scaled using and coad¬ 
ded, weighting each spectra by the rms noise. We rebin 
the spectra by a factor of 10, to achieve «14 kms“^ ve- 

^ http://www.iram.fr/IRAMFR/GILDAS 


locity bins. Typical rms per 14 kms ^ channel is 1 mK. 
These observations are summarized in Table |3l 

2 . 4 . Galaxy Properties from the SDSS 

We use a variety of data products from the SDSS to 
study properties of the post-starburst sample, including 
emission line fluxes, stellar masses, SFRs, and BPT clas- 
sifi cations from the MP A-JHU group catalogs (described 
in lAihara et al.l l2011h . We use Petrosian (jPetrosianI 
Il976ll optical sizes measured in the r band from the SDSS 
photometric catalogs, and redshifts from DR7. 

We use the stellar masses c alculated from the SPSS 
spectra (method described in iKauffmann et al.l l2003at 
iSalim et al.l l2007ll and included in the MPA-JHU data 
products. Because the star formation histories of post- 
starbursts may not be well represented by the templates 
assumed in the spectral fitting, we estimate the sys¬ 
tematic error by comparing stellar masses from several 
different algorithms run on SDSS data. We compare 
the stellar masses f rom the MPA-JHU data products to 
those calculated bylGhen et al.l (1201211. who rise both the 
iBruzual fc Charloti (|2003D and lMaraston fc StrombackI 
( 201111 SPS models. All three stellar mass calculations 
use a Kroupa initial mass function (IMF). The system¬ 
atic error from this method slightly exceeds the formal 
errors on the MPA-JHU measurements and is typically 
~ 30%. 

2 . 5 . Star Formation Rate Upper Limits 

We use two different methods to calculate SFRs for the 
post-starburst sample, one employing the Ha luminosity 
and the other the iJ„(4000) break. Both are contami¬ 
nated by other effects (principally LINER and A-stellar 
emission, see below), and serve as upper limits on the 
actual current SFR. 

Using the emission line f luxes from the MPA-JHU 
dataset (jAihara et al.l 1201 It) . we calculate SFR lim- 
its from Ha l umino sities using the relation from 
iKennicutt et al.l (|1994H . We use the Balmer decrement of 
Ha/H/3 to calculate dust extinction, assuming the stan¬ 
dard case B recombination at T= 10^ K and an intrinsic 
value of 2.86. We use the reddening curve of lO’Donnelll 
(|1994[) . For the cases where the H/3 line flux is uncer¬ 
tain, we use the mean value of E{B — V) of the other 
post-starburst galaxies. The mean attenuation is then 
Ay = 0.92 mag, or Ana = 0.77 mag. 

A complicating factor in determining the SFRs from 
Ha for the post-starburst sample is the high incidence of 
LINER spectra. A BPT diagram for the post-starburst 
sample is shown in Figure |TJ Two galaxies lie in the tran¬ 
sition region, and the rest are categorized as LINERs. 
Although the source of the LINER may not be an AGN 
(LINER emission is commonly seen in late st age mergers; 
Rich et all 120111 and from post-AGB stars; iSingh et"^ 
201311 . processes in addition to star formation will con¬ 
tribute to nebular line fluxes here, making the derived 
SFRs upper limits. 

The MPA-JHU group use the iJ„(4000) break as a 
less precise, but less contaminated way to estimate SFRs 
when galaxies do not l ie in the star-forming se quence 
on the BPT diagram (jBrinchmann et id] 1200411 . The 
iJ„(4000) break is a measure of the specific SFR (sSFR), 
and is calculated from regions of the rest-frame spec¬ 
tra bracketing the strong “break” observed near 4000A. 
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D„(4000) is not exp ected to be influenced by t he presence 
of a Type II AGN (|KaufFniann et al']l2003bD . The con¬ 
version between £)„(4000) and sSFR is calibrated from 
those galaxies in the SDSS categorized as star-forming. 
The scatter in this relation is large, and the error bars 
we show on the Z3„(4000) SFRs (derived using the MPA- 
JHU stellar masses) reflect the low precision of this cali¬ 
bration. 

The problem with using the Z3„(4000)-based SFRs in 
post-starburst galaxies is its sensitivity to the bright A- 
stellar populations produced in the recent burst. During 
ongoing star formation, the 4000A break is minimal, so 
D„(4000) is low. In passive galaxies, i4„(4000) is large. 
However, the timescale over which D„ (4000) is affected 
by a s trong burst (~1 Gyr, see e.g., iKauffmann et al.l 
l2003all is larger than the post-burst ages of our sample, 
so i4„(4000) here will reflect both previous and current 
rates of star formation. D„(4000) will be lower (more 
like star-forming galaxies) in post-starbursts than ex¬ 
pected given their instantaneous SFRs. Lower values of 
D„(4000) correspond to higher SFRs, so D„(4000) will 
overestimate the current SFR due to the recent burst. 
We use D„(4000)-based SFRs as upper limits on the cur¬ 
rent SFR. 

To use the SDSS fiber spectra to calculate global SFRs, 
we must account for any star formation outsid e the 3" 
fiber aperture. Like iBrinchmann et ^ (I2004f) . we see 
a trend of increasing fiber-based SFR per total stel¬ 
lar mass with redshift, after breaking up our complete 
SDSS post-starburst sample (1207 galaxies) into stellar 
mass bins. Thus, we expect some contribution to the 
SFR from outside of the fiber, so we require an aper¬ 
ture correction 0. We apply the aperture correction 


used in the MPA-JHU SFRs ( 

Brinchmann et al.l 12004 

iSalim et al.ll2007l lAihara et al.l 

201 ih. which is based on 


galaxy photometry outside the fiber. Although this aper¬ 
ture correction is calibrated on star-forming galaxies, it 
successfully removes the trend of sSFR with redshift for 
our complete post-starburst sample. While our use of 
this correction assumes that it also applies to our post- 
starbursts, the corrected SFRs remain likely upper limits 
as post-starbursts tend to have more positive color gra- 
dients (relatively bluer cores) than star-forming galaxies 
(lYang et al.llMlH [200l . 

One case where Ha and D„(4000) would not pro¬ 
vide upper limits on the SFRs is if we have significantly 
underestimated the dust extinction in post-starbursts. 
Radio continuum emission at 1.4 GH z is often use d 
as an “extinction-free” SFR indicator (|Gondonl Il992h . 
We search the FIRST ( Faint Images of th e Radio Sky 
at Twenty centimeters, iBecker e t al.l 1199 511 an d NVSS 
(NRAO VLA Sky Survey. IGondon et al.l I199§1 1.4GHz 
surveys for matches within 10" of each galaxy in our 
sample. We find 6 detections in the FIRST survey (S06, 
HOI, H03, H07, H08, and H09). The galaxies H07, H08, 
and H09 are als o detecte d in th e NVSS. Using the con¬ 
version found in iGondonI (jl992tl . the SFRs suggested by 
these detections are higher than the Ha SFRs. If we 
were to accept that the standard 1.4GHz - SFR relation 
is valid for the post-starburst sample, it would require 
these galaxies to have up to 4.7 magnitudes of addi- 

® We consider the case where star formation is limited to the 
fiber aperture in E21 see Figure Ish. 


tional extinction on top of the ^ 1 magnitude already 
accounted for using the Bahner decrement. While dust 
extinctions of 5-6 magnitudes are not unheard of, es¬ 
pecially for starbursts, the dust extinction is consistent 
with that implied by Balmer decremen t in those cases 
(jGhoi et al.ll2006t iKennictitt et al.ll2009| l. The huge dif¬ 
ference between the extinction derived from the Balmer 
decrement and implied by the 1.4Ghz-SFR relation is 
unprecedented and suggests a problem with the SFRs 
derived from the 1.4GHz data for our post-starbursts. 

The LINER and recent starburst in these galax¬ 
ies complicates the standard 1.4GHz SFR conversion. 
Galaxies with LINER spectra have enhanced 1.4GHz 
fluxes when compared to other measure s o f their SFRs 
( de Vries et al.l l2007t iMoric et al.l 1201011 . iMoric et al.l 
( 2010(1 find that 90% of the 1.4GHz flux can come from 
the LINER, not from star formation, and that the scat¬ 
ter in the 1.4GHz-SFR relation for LlNERs is large, of 
order 2 dex. For our sample, the 1.4GHz-based SFRs 
and limits scatter evenly about the Ha-based SFRs af¬ 
ter the radio SFRs are reduced by the expected factor 
of 10. While the LINER will also contribute to the Ha 
flux, its contribution is ty pically < 40%, with less scatter 
(|Brinchmann et ^120041 1 , implying that Ha is more reli¬ 
able than 1.4GHz as a SFR upper limit. Additionally, the 
recent large starburst may boost the amount of 1.4GHz 
flux on time scales overlapping wi th the post-burst ages of 
our s ample (jBressan et al.l[200^ . As discussed above in 
1 12.11 the galaxies marked “S” were selected with a cut on 
the [0111] equivalent widths, intended to exclude galax¬ 
ies with strong AGN activity from the sample. This cut 
was not applied to the selection of galaxies marked “H,” 
and the higher incidence of 1.4GHz detections in the “H” 
sample may be tied to a higher incidence of AGN. 

Both the TIR lum inosity (total I R, from 8-1000 fxm ; 
iHavward et al.l[2?)l^ and 24 ^m flux (jUtomo et al.ll20l3l 
are strongly affected by dust heating by the substantial 
A stellar population in post-starbursts, so we do not con¬ 
sider these SFR indicators here. Ongoing analysis of our 
sample observed in PAH emission and high ionization 
species (Smercina et ah, in prep) will provide further 
constraints on any current SFR. 

In the following analysis, we use the Ha-derived SFR as 
an upper limit, as well as showing the effect of assuming 
the Zl„(4000)-based SFRs. 

3. RESULTS 

3.1. Detection of Molecular Gas 

We detect molecular gas at > 3cr in 17 of the 32 galax¬ 
ies observed, using IRAM 30m measurements of the GO 
(1-0) line. If we increase our detection threshold to 4cr, 
we detect 14 galaxies, and at > 5a, we detect 11 galaxies. 
To calculate the integrated CO line intensity Ico, we fit 
a Gaussian profile to each line, allowing the center veloc¬ 
ity to differ from the optical velocity up to 200 kms“^. 
We use the Gaussian width Cgauss to choose integration 
limits of iSCTgauss- Although many of the line shapes 
are not exactly Gaussian, this method allows us to es¬ 
timate appropriate velocity intervals for integration in a 
systematic way. FWHMs given by these fits are listed in 
Tables El and 131 If the signal to noise ratio for Ugauss is 
< 3, we use the interval [-260, 260] kms“^ These veloc¬ 
ity limits were chosen to be the median of those of the 
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log([N II] A 6583)/(Ha A 6563) 

Figure 1. BPT JBaldwin et aljfl98U) diagram for post-starburst 
sampl e, measured from SPS S spectra. Galaxies from the SDSS 
DR7 IIAbazaiian et ani20Q9l) with well-measured lines are shown 
as a shaded backgrou nd. The lines separa ting star-forming and 
AGN-li ke activity from fl^wlev et al.l II2001J) and IKauffmann et al.l 
(l20Q3bl) are shown as dotted and solid lines, respectively. The line 
at $ = 25 degrees separates Seyferts from LINERs. The post- 
starburst sample is plotted as individual points, for galaxies with 
all lines detected at > 3cr, with a characteristic errorbar shown 
in the bottom right. Most of the post-starburst sample, except 
S12, is solidly in the LINER category. The presence of LINERs 
complicates our calculation of the current SFR, as the nebular 
emission lines will be contaminated. 


well-fit sample, and are centered around the optical ve¬ 
locity. The velocity intervals fit from the CO (1-0) data 
are used for the CO (2-1) data, though we note that fit¬ 
ting the CO (2-1) data separately does not change our 
results by > Icr. We calculate the error in the integrated 
CO line intensity as 

aj = {Avf cr^ Ni (1 -h ^), (1) 


where Av is the channel velocity width, a is the channel 
rms noise, Ni is the number of channels used to integrate 
over the line, and Nb is the number of channels used to 
fit the baseline. We also take into account an estimated 
flux calibration error of 10%. We calculate upper l imits 
on I CO as < 3(77. Following ISolomon et al.l f|1997ll . the 
CO line luminosity L'qq (in Kkms“^pc^) is 

L'co = 23.5 n^^b Dl Ico (1 + ^)“" , (2) 

where Ico = / Tmb dV is the integrated line intensity 
(in Kkms“^) as described above, z is the SDSS redshift, 
and Dl is the luminosity distance (in Mpc). is the 
solid angle of the source convolved with the beam, 


n 


s*b — 


4 In 2 


(3) 


where 9s and 9b are the half power beam widths of the 
source and beam, resp ectiv ely. Because the CO emit¬ 
ting size estimates (see M.2D are not available for all the 
sources, we adopt a simple approximation such that the 
beam is much larger than the source, so Ils*b ~ Ufa- Note 
that depending on the actual size estimates in M.21 we 
could be underestimating L'co by ^ 1.1 — 2.2x, with a 
median of 1.4, but this does not affect our conclusions 
throughout the paper. 


The molecular gas mass can be calculated from 
by assuming a conversion factor aco y as 

M(H2) = acoL'co- (4) 

For now, we assume an aco comparable to that in 
Galactic molecular clouds and the Local Group (aside 
from the SMC; see recent reviews bvIBolatto et al.l 120131 : 
iCarilli fc WalteilF2013HCasev et al.ll2014l l: aon = 4 Mq 
(Kkms“^ pc^)“^ (units omitted hereafter). This choice 
is examined below and in ii4.3l 

Molecular gas masses for the post-starburst sample 
span a broad range, from 3.4 x 10® to 6.9 x 10^ Mq, with 
a mean value of 3.0 x 10 ®Mq among the detected sample. 
We measure upper limits for the remaining 15 galaxies, 
with 3(7 limits ranging from 4.6 x 10^ to 5.2 x 10®M©. 
Molecular gas masses and upper limits are listed in Ta¬ 
ble [21 Optical postage stamps of the galaxies with and 
without molecular gas detections are shown in Figures [T] 
andini respectively. 

Next, we compare the molecular gas masses measured 
here to those from surveys of other galaxy types. CO 
(1-0) measurements have been compiled for the At las-3D 
sample of earl y type galaxies (Youiig et al.ll201lD . The 
COLD GASS ([Saintonge et al.l 12011 ) sample is a stellar 
mass-limited sample of galaxies, selected from the SDSS 
independent of galaxy type. We divide the COLD GASS 
sample up by galaxy type assigned by the SDSS based 
on the optical spectra (galspec bptclass). For now, we 
only use galaxies classified as star-forming or low signal- 
to-noise star-forming. We assume aco = 4 to calculate 
molecular gas mass for the early type and COLD GASS 
star-forming samples. 

We compare the total molecular gas masses of these 
samples in Figure lU seeing significant overlap between 
the star-forming and post-starburst samples. This over¬ 
lap is surprising, because of the lack of equivalent levels 
of current star formation in the post-starburst sample. 
The lower mass end, as well as the upper limits, of the 
post-starburst sample are consistent with the early type 
sample. 

In addition to comparing M(H 2 ), we also compare 
molecular gas fraction /gas = M(H 2 )/M* normalized by 
stellar mass M*. We us e calculated from the SDSS 
spectra, as discussed in H2.41 for both the post-starburst 
and COLD GASS samples. We calculate stellar masses 
for the early type galaxies in the same way as Atlas -3D, 
using K-band measurements (jCapDellari et al.ll2011h . 

As with the M(H 2 ) comparison, we see considerable 
overlap in /gas between the post-starburst and star¬ 
forming samples. These comparisons are shown in the 
right-hand panels of Figure [2] The molecular gas frac¬ 
tions for the post-starburst sample are primarily above 
those of the early type sample, while some of the upper 
limits are more consistent with early types. 

The CO to H 2 conversion factor (aco) is a known 
source of uncertaint y in observations o f molecular gas 
(see recent review bv IBolatto et al.ll2013l l. Traditionally, 
a bimodel model has been used, with normal star-forming 
galaxies assigned a Milky Way-like value of aco ~ 4, and 
ULIRGs or starbursting galaxies assigned aco 0.8. 
This approach was initially motivated by the fact that a 
high aco applied to ULIRGs produced gas masses higher 
than the dynamical masses. In Figure |3l we plot M(H 2 ) 
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Figure 2. Left: Histograms of derived molecular gas mas ses M(H 2 ) for a vari ety of galaxy types: early types (top, from Atlas-3D 
I Young et al.|[20T^^ . star-forming (bottom, from COLD GASS: [Saintonge et al.ll201111 . and our post-starburst sample (middle), aco = 4M0 
(Kkms“^ pc^)“^ is assumed for all samples. Bin size represents the mean error in the post-starburst sample, excluding systematic error 
from uncertainties in oco- A histogram of 3cr upper limits is overplotted tor non-detections. Right: Histograms of molecular gas normalized 
to stellar mass (/gas) for the same samples. For both M(H 2 ) and /gas, we see considerable overlap between the post-starburst sample and 
star-forming samples, which is surprising given the difference in SFRs. The lower end and upper limits of the post-starburst sample are 
consistent with AL(H 2 ) and /gas measured for the early type sample. As seen in Figure [S] overlap persists even if a ULIRG-type value of 
OCO is assumed for the post-starburst sample. 


and /gas for the post-starbursts and comparison galax¬ 
ies for different aco assumptions. We also compare to 
the sa mple of LIRGs and ULIRGs from iGao fc Solomonl 
(|2004D . Even if a low, ULIRG-like value of aco = 0.8 is 
used for the post-starburst sample, we still see significant 
overlap with the star-forming sample (at aco = 4), and 
even some overlap with the LIRG and ULIRG sample 
(at aco = 0.8). We expect aco = 0.8 and aco = 4 to 
span the range of possible values of aco in post-starburst 
galaxies, because recently ended starbursts may reflect 
ISM physical conditions between ULIRGs and quiescent 
disk galaxies. However, the appropriate value of aco 
for post-starbursts remains largely unconstrained. We 
discuss the effects of this uncertainty on our results in 

SH 

3.2. High Molecular Gas Mass for Given SFR 

We compare the molecular gas masses measured here 
to upper limits on SFR derived from Ha and i7„(4000) 
in Figure [B] We use several comparison datasets: 
star-forming galaxies select ed from the GOLD GASS 
(|Saintonge et al.ll201lL[201^ . the star-forming and star- 


burst sample from IGao fc Solomonl (120041). and early 
type g a laxies (with Ha-|-P AH SFRs) from lYoung et al.l 
(|20Ili) : iDavis et'all 1)20141) . The post-starburst sample 
lies at higher M(H 2 ) for their SFRs than early-types, 
star-forming galaxies, and (U)LIRGs. The median Ha- 
derived SFR upper-limit for the post-starburst sample is 
~ lOx lower than the median SFR for the COLD GASS 
star-forming sample across the same A/(H 2 ) range. This 
offset persists for the median of the D„(4000)-derived 
SFR upper limits, which is 20x lower than expected 
given M(H 2 ). 

Is it possible that our post-starburst selection crite¬ 
ria have generated the observed offset in SFR? Because 
we selected the post-starburst sample to have low Ha 
equivalent widths, it may include star-forming galaxies 
whose EW(Ha) measurement errors have scattered them 
low. This does not appear to be the case, as star-forming 
galaxies (classified using bptclass=l) that pass our H^ 
absorption cut typically have EW(Ha)^ SA, the limit 
for our post-starburst sample. A Monte-Carlo analysis 
predicts a 0.004% contamination rate, and even if the 
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Figure 3. Left: Histograms of derived molecular gas masses 
M(H 2 ) for a v ariety of galaxy types: star-forming (bottom, from 
COLD G ASS: [ Sainton^e et al.|| 2011l) . LIRGs and ULIRGs (bot¬ 
tom, from IGao omoiJT2004 ). and our post-starburst sample 

(top), aco = 4Mq (Kkmpc^)“^ is assumed where data are 
plotted as solid histograms, and oco = O.SMq (K kmpc^)“^ 
where histograms are dashed lines. Bin size represents the mean 
error in the post-starburst sample. Right: Histograms of molecular 
ga s normalized to stellar mass (/gas) for the same samples (except 
for IGao fc Solomonll200^ . Even if a low, ULIRG-like value of oco 
is used for the post-starburst sample, we still see significant over¬ 
lap with the star-forming sample, and even some overlap with the 
LIRG and ULIRG samples. 

Ha equivalent widths had systematic errors 3x as large 
as their measurement errors, we still expect <C 1 con¬ 
taminant in the complete SDSS post-starburst sample of 
1207 galaxies. Thus, the KS offset does not arise from 
contamination from normal star-forming galaxies. 

The large molecular gas reservoirs in post-starburst 
galaxies are inconsistent with their SFRs when compared 
to a broad sample of galaxy types. Thus, the cessation of 
star formation after the starbust cannot be due to a lack 
of gas in the nearly half of our sample with detected CO. 
The question remains: why are these galaxies no longer 
forming stars at significant rates? One possibility is that 
the molecular gas is spread out over a larger area, drop¬ 
ping its surface density to be consistent with SFR den¬ 
sity o n the Kennicutt-Schmidt (KS) relation (iKennicutd 
119981 hereafter K98). We examine the KS relation below. 

3.3. Offset from the Kennicutt-Schmidt Relation 

While there are clear trends of SFR with molecular gas 
mass, tighter correlations exist when comparing the sur¬ 
face densities of these quantities for normal star-forming 
and starburst galaxies. 

We determine the molecular gas surface density Ehj 
and SFR surface density Esfr for the post-starburst 
sample, using the SDSS r-band Petrosian 90% radius to 
calculate the area as 

Eh 2 = acoEco/(^-^9o); ^SFK = SFR/ {ttRIq). (5) 

We place the post-starburst sample on a KS plot in Fig¬ 
ure [71 Other local galaxies are shown for comparison, in¬ 
cluding both normal star-forming galaxies and starbursts 
from the canonical K98 sample. For now, we apply the 
same value of aco = 4 to the post-starburst sample, 
and the entirety of the K98 sample. Many of the post- 
starburst galaxies lie below the relation defined by the 


other galaxies: the Ha-derived SFR limits exclude con¬ 
sistency with the relation for all but 4-5 galaxies. The 
D„(4000) based SFR limits also lie mostly below the re¬ 
lation. 

The median locus of the 17 post-starburst galaxies lies 
gX lower than the n = 1.4 power law fit to the K98 
galaxies. We perform a Monte Carlo analysis to assess 
the significance of this result by choosing random sets 
of 17 galaxies from the K98 disk sample, finding a 5a 
significant offset for the post-starbur st locus. Thi s offset 
is more extreme than that found by iDavis et al.l ()2014r i 
for their sample of gas-rich early type galaxies. The re¬ 
lationship between the two datasets is unclear. 

We see no obvious differences between the properties 
of our galaxies that are roughly consistent and most dis¬ 
crepant with the KS relation. It is not clear if the post- 
starbursts are single population of galaxies, or several 
different families. 

The optical size used by K98 to calculate the surface 
densities is the isophotal radius, where the B-band sur¬ 
face brightness drops to 25 mag arcmin“^ and which is 
comparable to the Ha emitting region for normal spi¬ 
ral galaxies (K98, although we test this assumption for 
our sample in >14.21) . This isophotal radius should be a 
good estimate of the size if the CO emission is coming 
from the same region as the optical light from star for¬ 
mation. Here, we use the Petrosian 90% radius for our 
post-starbursts, because the isophotal sizes in SDSS are 
not considered reliable and are not included in the pho¬ 
tometric catalogs after DR8. However, the significant 
offset from the K98 sample remains if we use the r-band 
isophotal radii instead. 

There are several other observational uncertainties 
that could affect our results, which we consider in the 
Section 4. We d iscuss the effect of our sample selection 
criteria in SH Like K98, we assume that the optical 
size of the galaxies is a good proxy for the spatial extent 
of both star formation and molecular gas. This assump¬ 
tion may not be valid for post-starburst galaxies. In >14.21 
we test the possibilities that the CO is distributed differ¬ 
ently from the optical light and that the star formation 
is distributed differently from both the optical light and 
most of the H 2 , as traced by CO. We also consider how 
our sampling of the CO region (aperture bias) might af¬ 
fect our results. In >14.31 we test for the possibility that 
the measured CO is not tracing the H 2 as we expect, 
resulting in a different aco- 

4. DISCUSSION 

4.1. Possible Sample Selection Biases 

Given the way we selected our CO targets, the sample 
observed here might not represent the gas properties of 
the overall post-starburst sample. To study any biases 
that may occur within our sample, we test whether the 
galaxies with CO (1-0) detections lie at the extremes of 
our selection criteria. 

The two parts of our sample (labeled “H” and 
“S”) were selected from the parent sample of SDSS 
post-starbursts using different criteria (more details in 
Smercina et ah, in prep). The “H” sample was selected 
based on post-starburst galaxies bright in the WISE 
12^m band. One might expect these galaxies to have 
more gas if the 12^m band is a proxy for hot dust con- 
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Figure 4. 60" by 60" SDSS postage stamps of the 17 post-starburst galaxies with CO (1-0) detected at > 3cr with the IRAM 30m. The 
size of the 3" SDSS fiber is overplotted in blue, and the size of th e IRAM 30m CO (1—0) 22" beam is overplotted in orange. Galaxies 
are ordered by decreasing M(H 2 ). Given the size estimates in N4.2I we could be underestimating L'qq by factors of 1.1 — 2.2x, with a 
median of 1.4, due to aperture effects. 
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Figure 5. Same as Figure[4] but for the galaxies not detected in CO (1-0). 
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Figure 6. Molecular gas mass vs. SFR for post-starburst (PSB) galaxies and comparisons. Whether Ha (left) or Dti( 4000) (right) SFR 
upper limits are used for the post-starburst sample, these galaxies fall s ystematically below the comparison galaxies from t he COLD GASS 
sampl e classed by the SDSS as star-forming IISainton^e_etal .[200. [201211 . star-forming, LIRG, and ULIRG galaxies from the lGao &: SolomonI 
(1200411 sample, and early type galaxies from IYoTing^^dTluOM^ . Both Dn(4000) and Ha are expected to overestimate the SFR in the 
post-starburst sample, so these galaxies may lie at even lower SFRs. All galaxies have been normalized to the same value of aco = 



Figure 7. Molecular gas surface density vs. SFR su rface density from Ha (left) and Z)n(4000) (right) SFR upper limits. The post-starburst 
sample is shown in black, with other local galaxies IIKennicutUll99^ shown in red. Red circles are normal star-forming galaxies and red 
squares indicate local starbursts. Both SFR indicators are upper limits to the true current SFRs for the post-starburst sample. Ha has a 
contribution from the LINER, and Dn(4000) has a contribution from the intermediate age stars produced in the burst. Our post-starburst 
galaxies are biased low on this relation. For the post-starburst sample, we assume the CO is distributed in the same way as the star 
formation regions, using the Petrosian 90% radius Rgo to calculate the surface densities, EgpR = SFR/tvRqq and Eh, =M{H2)/nRlo. 
Here, M(H 2 ) includes all the CO detected in the IRAM 30m beam. We use the same value of aco = 4 for all galaxies. This plot is the 
most appropria te co mparison to the K98 data set, although we test our assumption that CO and star formation are distributed like the 
optical light in E2] and that aco = 4 in E31 
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tent and their dust traces their gas. However, we see 
no mean offset in 12/im luminosity between the galaxies 
detected and not detected in CO. 

Both samples were selected to have a variety of times 
elapsed since the starburst (post-burst age). If the molec¬ 
ular gas is depleted over time, younger post-starbursts 
may be easier to detect in CO. While there is a shift to¬ 
wards detections with younger post-starburst ages, it is 
not statistically significant. The relation between molec¬ 
ular gas content and age since the burst is not straight¬ 
forward, and is heavily dependent on the pre-burst gas 
mass of the galaxy and on the mechanics of the burst 
itself. 

Because no statistically significant boost in CO detec¬ 
tions occurs with either younger post-burst age or higher 
12 ^m luminosity, the molecular gas properties of our 
sample here are not significantly biased by the selec¬ 
tion criteria. Therefore, our sample is likely to have a 
molecular gas detection rate representative of the overall 
post-starburst population. 

4.2. Effect of Spatial Distributions of Gas and Stars 

In ^3.31 we made the assumptions that both the CO 
and current star formation were well-traced by the op¬ 
tical light, in order to calculate their surface densi¬ 
ties. ldn^_ 2 S_a^ood_assum 2 tion_for star- forHung_gaInr^ 
ies (iRegan et al.l 120011 : iLerov 611111120081 iSchruba et al.l 
1201 ih . but may not apply to our sample, especially if the 
reason for the end of the starburst is a disruption of the 
gas. Additionally, the spatial extent of any residual star 
formation is unknown and may not overlap with the op¬ 
tical light, which is dominated (~ 60 — 90%) by the A 
stars formed in the recent burst. Without resolved ob¬ 
servations, we are limited in how accurately we can know 
the distributions of gas and current star formation. If the 
CO emission or any currently star forming regions have 
sizes different than the optical size, it might be possi¬ 
ble to resolve the observed offset from the KS relation. 
First, we test these assumptions for the post-starburst 
sample by estimating the CO emitting size with a model 
for the CO emission. We continue here to assume that 
CO traces H 2 well, with a conversion factor of aco = 4. 
Second, we use our estimate of the CO emitting size to 
compare the scaled amount of CO near the center of the 
galaxy to the SDSS fiber-based SFRs. 

By comparing the CO emission in two differently sized 
beams, we can roughly constrain the CO emitting size by 
assuming a Gaussian model for the shape of the emitting 
region. We estimate the source size for each galaxy by 
combining the CO (2-1) line measurement s from IRAM 
30m a nd SMT 10m, using the method from iLavezzi et ahl 
(fT^ . IcOi the integrated line intensity, is related to 
the surface brightness, so it should scale with the con¬ 
volved size of the beam and the source as 

ISMT _ + ^b.IRAM 

IlRAM 4 " ^b,SMT 

for the same line, assuming a Gaussian distribution of 
the CO emission, where 9i, are the different beam sizes 
for IRAM 30m and SMT, Os^smt = 1.2X/D « 33 arcsec, 
and Ob.IRAM = 1.166A/I? « 11 arcsec, with the different 
coefficients due to the taper of each telescope. The source 


size Og is then given by, 

I I RAM (^bJRAM ~ ^SMT (^b,SMT 
I SMT — IlRAM 

for a Gaussian source 

We estimate source sizes for the 6 galaxies with IRAM 
30m CO (1-0), CO (2-1) and SMT CO (2-1) detections. 
Because the measurement errors propagate non-linearly, 
we use a Monte Carlo technique to estimate the formal 
error on the size estimates, excluding any systematic er¬ 
ror from the Gaussian model assumption. Sizes are listed 
in Table 131 and plotted over SDSS postage stamp images 
in Figure [ 8 ] 

Previous studies of star-forming galaxies have found 
comparable exponential scale lengths of the CO emission 
and optical emission u s ing resolved data (IReean et al.l 
120011 : ILerov 611111120081 ISchruba et al.l 1201111 . However, 
the CO to optical sizes are on average ^ 2x larger for 
our post-starburst sample, using either the optical half- 
light radii from an exponential fit or the Petrosian half- 
light radii. Because of the centrally peaked distribution 
of A stars in the post-starburst sample, this difference 
is due to the concentrated optical light, rather than ex¬ 
tended CO emissio n. Concentration s measured for post- 
starburst galaxies (|Yang et al.ll2008ll are high, consistent 
with the half light radii being smaller than for galaxies 
with an exponential profile. 

To test the effect of assuming a Gaussian model, we 
perform a similar calculation for each galaxy, but for 
a uniform disk emitting region instead of a Gaussian. 
The half-light sizes for each source model are consistent 
within the measurement errors except for EAS06, where 
the uniform disk model predicts a size smaller than the 
Gaussian prediction (20" vs. 25" ). 

The systematic errors associated with this method may 
be significant, especially as several of the estimated sizes 
are larger than the CO (2-1) IRAM 30-m beamsize. 
Thus, we use these sizes only to roughly estimate aper¬ 
ture bias in the CO observations, and to test whether use 
of this alternate size measure can eliminate the offset in 
the observed KS relation. The CO size estimates play no 
role in our main conclusions. 

We use the CO emitting size estimates to model the ef¬ 
fect of aperture bias in our CO (1-0) observations taken 
with the IRAM 30-m k,22'' beam. Given that the CO- 
emitting region is not a point source, the CO line lumi¬ 
nosity calculated in §3.1 is likely an underestimate by a 
factor of (0g -I- 0l)l9"l, where 9s and 9b are the source 
and beam sizes, respectively. If we assume the mean CO 
source size of 9.2 kpc from our estimates above, then 
the corrections needed for the CO line luminosity, thus 
also M(H 2 ) are ~l.l-2.2x. The mean correction is an 
increase of 1.4 x from the original calculation, indicat¬ 
ing that the molecular gas masses reported here may be 
conservative, with even higher amounts of molecular gas 
present. 

Using the models of the CO emitting region that we 
construct from the size estimates and Gaussian assump¬ 
tion, we test the possibility that some of the observed 
CO does not participate in star formation. In il3.3l and 
in Figure [71 we assumed that all of the observed CO 
lay within the optical size R 90 , which is always smaller 
than the IRAM 30m CO (1-0) beamsize. Our CO size 

























12 



Figure 8. Estimated Gaussian half-light sizes of CO emitting 
regions (white dotted circles), with horizontal lines representing 
the Monte Carlo estimated error. FWHM beam sizes are shown 
for comparison: the inner yellow circle represents the IRAM 30m 
CO (2-1) Rill” beamsize, and the outer orange circle is the SMT 
CO (2-1) Ri33” beamsize. Optical images are from SDSS, and are 
60” X 60” . These CO sizes are ~ 2x the optical half-light sizes 
(not shown), due to the concentrated optical profile from the stars 
produced during the burst. 

estimates, however, indicate that CO may extend be¬ 
yond this aperture. We now estimate Sh 2 within RgQ 
by rescaling M^Hg) using the Gaussian model to include 
only the gas mass within Rgg (Figure IH^). While still as¬ 
suming that the current star formation is distributed like 
the optical light, this method allows us to compare the 
surface densities within the same aperture. In il3.1l we 
assumed that the source sizes were much smaller than the 
beam sizes to calculate ^( 7 o(i_o)i but we now calculate 
^ls*b explicitly in Equation 1. Despite the decrease in 
M[H 2 ) (and the post-starburst sample is still sig¬ 

nificantly (> dtr) offset from the K98 galaxies. Allowing 
for the optical light and CO to be distributed differently 
does not change our results. 

Next, we study the effect of allowing the spatial extent 
of any current star formation to differ from both t he CO 
and optical light distributions. iSwinbank et al.l (j2012f ) 
use resolved IFU observations of post-starburst galaxies, 
and find the nebular emission lines [O II] and [O III] to 
be spatially offset from the A star population in some 
cases, although these lines are particularly contaminated 
by LINER emission. If the current star formation is even 
more concentrated than the optical light, this could drive 


the observed offset of the post-starbursts on the KS rela¬ 
tion. If this is the case, our assumptions about the aper¬ 
ture correction in the SFRs in 112.51 are not valid, and the 
SFRs are even more of an over-estimate. We consider 
now the case where the star formation is restricted to 
the optical fiber. 

Using the CO sizes estimates and Gaussian model, we 
rescale M{H 2 ) to that within the 3" SDSS optical fiber, 
again with a Monte-Carlo method to estimate the errors. 
We then estimate Ehs inside the hber aperture by scal¬ 
ing the CO luminosity, assuming a Gaussian profile with 
the GO size estimate. To calculate Esfr , we use only 
the Ha flux from the fiber, not the aperture-corrected 
flux used elsewhere throughout this paper. We plot the 
resulting surface densities within the fiber apertures in 
Figure [5 Jd. The post-starburst galaxies are still signifi¬ 
cantly offset from the K98 comparison galaxies. Thus, 
even if we allow for the possibility that some of the GO 
in the outer regions of the galaxy does not participate 
in star formation, or that the current star formation is 
more compact than the optical light, the offset from the 
K98 galaxies remains. 

Even if the CO in the post-starburst sample and the 
K98 sample is distributed in the same way with respect to 
the optical light, aperture bias in the datasets could gen¬ 
erate an apparent offset. We do not see evidence of severe 
aperture bias in the post-starburst sample. Although 
we do observe higher molecular gas masses at higher 
redshifts, this is due to a combination of our decreased 
sensitivity and the higher stellar mass SDSS-selection at 
higher redshifts. We see no statistically significant trends 
with redshift of /gas, or in the offset from the n = 1.4 
power law KS relation from K98. If aperture bias only 
affected the post-starburst sample, it would result in an 
under-estimation of the offset. 

An offset could be generated, however, if the K98 sam¬ 
ple was not measured out to the same physical radii as in 
the post-starburst sample. The CO isophotal sizes in the 
K98 star-forming sample are at most 60% of the optical 
isophotal sizes, and often much smaller. Galaxies with 
CO measured out to at least half the isophotal radii are 
the only ones included in the K98 sample. Assuming a 
worst-case, uniform disk distribution, the actual CO flux 
could be up to 70% higher than measured, which is not 
enough to resolve the observed offset of ~ 400%. 

The offset observed in Figured that the post-starburst 
sample lies at lower SFR for a given M(H 2 ) than other 
galaxies, and the similar robust offset in the KS relation 
(Figures [7] and ini), suggest that I) the CO does not trace 
the H 2 as in our comparison galaxies or 2) the H 2 is not 
turning into stars in the same manner as our compari¬ 
son galaxies (i.e., the star formation efficiency (SFE) is 
lower or the IMF is bottom-heavy). We explore I) in 
Section and 2) in Section For now we note that 
suppressing the SFE by allowing CO to extend beyond 
any current star formation region does not resolve the 
KS offset. Ultimately, interferometric CO observations 
will be required to test the spatial distribution of CO. 

4.3. Effect of aco Choice 

As discussed in §3, the CO to H 2 conversion factor 
(aco) is a known source of uncertainty in measuring 
molecular gas masses from CO observations. Here, we 
assume that H 2 traces the current star formation as ex- 
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Figure 9. Alternative KS plots (red comparison points same as in Figure[7]|, accounting for possible spatial differences between molecular 
gas, star formation, and optical light. Left: Sspu (from Ho) and are both calculated within the optical radius Rgo for the post- 

starburst sample. We use the CO size estimate from H4.2l and assume a Gaussian distribution to rescale M(H 2 ) to that within Rqq instead 
of within the larger IRAM 30m 22” beam as in Figure [T] Despite this rescaling, which allows for the possibility that the optical light and 
CO are distributed differently, we still observe an offset of the post-starbust galaxies from the K98 sample. Right: Sgpjr and Sjjj are 
both calculated within the 3” SDSS fiber. We rescale M(H 2 ) to within this radius. SFR is calculated using only the Ha flux from the 
fiber, without rescaling. This is a test of whether the star formation is distributed differently than either the optical light or CO, and the 
continued offset shows that this possibility is unlikely unless the Gaussian assumption for the CO distribution is poor, or the star formation 
is even more concentrated within the fiber size. 


pected for other galaxies, and explore variations in how 
CO traces the H 2 . 

So far, in plotting the KS relation we have assumed a 
single value of ctco = 4 for all samples. Now assuming 
a bimodel aco model instead, we apply a ULIRG-like 
value of oco = 0.8 to the K98 starbursts and to our 
post-starburst sample, and leave the K98 star-forming 
galaxies with aco = 4. We obtain the results in Figure 
[inii. This low value of aco applied to the post-starburst 
sample can remove their observed offset from the modi¬ 
fied KS relation. 

We can understand why starbursting galaxies may re¬ 
quire a lower value of aco using the following toy model, 
aco is proportional to the column density of molecular 
gas N{H 2 ) over the CO line intensity, as 


aco C!C 


N{H2) 

Ico 


N{H2) 
T X a 


( 8 ) 


(jNaravanan et al1l2011t) . In a merger, the column density 
is increased. However, the line intensity goes as the tem¬ 
perature T times the velocity dispersion cr, which both 
increase during the merger. In total, these factors re¬ 
sult in a lower value of aco- After the merger, the gas 
kinetic temperature may decline to match conditions in 
early-type galaxies, but it is not clear what this simple 
model predicts for the post-starburst sample. 

We cannot assume the post-starbursts will have a simi¬ 
lar value of aco as ULIRGs or starbursts simply because 
they are the likely progenitors. We expect the physical 
state of the gas to have changed significantly in the 0.3- 
I Gyr since the burst, as the dynamical timescales for 


ULIRGs are of order 10® — 10^ ® yr (|Genzel et al.ll20i^ . 
However, if the state of the gas after the merger changes 
in such a way as to keep the gas heated, but as to lower 
the column density of molecular gas, aco may have a 
low value. Additionally, if an AGN heated the bulk of 
the molecular gas, the GO brightness temperature could 
increase for a given H 2 mass, lowering aco- 

If much of the GO emission comes from outside of 
GMGs, the GO (1-0) linewidth could be strongly af¬ 
fected by the gravitational potenti al in the galaxy, in¬ 
stead of just by its own turbulence (jPownes fc SolomonI 
Il998f) . The fact that starbursts have more diffuse gas, 
and less gas bound in molecular clouds, has been used 
as justification for their low values of aco- We can esti¬ 
mate the influence of the stellar pot e ntial o n aco using 
the prescription from iBolatto et aP ()2013[ ). They sug¬ 
gest that aco scales down from the Galactic value as 
acolaco.MW = (Stotai/100 Mqpc”^)”® ®, where Stotai 
is the combined surface density of stars and gas. Adding 
the stellar mass to the molecular gas mass, and calcu¬ 
lating the surface density within the optical radius i? 9 o, 
our post-starbursts have total surface densities of 130-460 
Mcpc”^, predicting aco = 1-9 — 3.5, scaled down from 
the Galactic value of aco = 4. This difference between 
the assumed aco for comparison star-forming galaxies 
and the post-starburst sample is not enough to resolve 
their observed offset. However, this relation is subject 
to scatter, and there may be si gnificant variation i n the 
influence of diffuse gas on aco (|Liszt fc Petvll201^ . 

The use of a bimod el aco is not necessarily physi- 
cal (although see, e.g., iTacconi et al.ll20()8l : iDaddi et al.l 
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[20n9h . because a variety of ISM conditions should exist, 
resulting in a continuum of aco values. In particular, 
because post-starburst galaxies may be at an intermedi¬ 
ate stage between being dynamically hot and more re¬ 
laxed, the appropriate aco for these systems may lie 
between the Galac tic average and th e typic al value used 
for ULIRGs. The iNaravanan et al.l (l2012f) formulation 
suggests that aco varies smoothly with galaxy physical 
properties, and may be parametrized in terms of the gas 
phase metallicity and the CO surface brig htness . Our 
galaxies do not have abnormal metallicities. iGotol (j2007f ) 
study the metallicities of post-starburst galaxies and do 
not find them to be anomalo us. The mass-metallicity 
relation (|Tremonti et al.ir2004D also predicts metallicities 
that are similar to the comparison galaxies, so we vary 
aco only with CO brightness. The result of applying 
this variable aco model to both the post-starburst and 
K98 samples can be seen in Figure ITOb . The variable 
aco model does not remove the offset. 

We are unable to rule out a low value of aco = 0.8 as a 
potential explanation for the offset from the KS relation. 
We can resolve this question observationally using higher 
Jup lines of CO to constrain the temperature and density 
of the gas, as well as with denser gas tracers such as 
HCN to probe denser regions of the gas, bypassing the 
uncertainties arising from any CO outside of GMCs. 

4.4. Implications for Galaxy Evolution 

Our findings that post-starburst galaxies can have 
large gas reservoirs, and that they are offset low from 
the KS relation, help to discriminate among the physical 
processes proposed to end the burst. Clearly, scenarios 
that require the molecular gas to be absent, such as the 
complete 1) expulsion or removal of the g as (e.g., in out¬ 
flows or some environmental mechanisms: iFeruglio et al.l 
120101: iRuoke fc Veilleu^ 120111 : iBoselli fc Cavazzil 120061 1. 
2) consumption of the gas in the burst, and 3) prevention 
of gas infall into the galaxy and of the subsequent for - 
mation of new GMCs 1 “starvation”: iLarson et anil98ClD . 
are now excluded, at least in the half of our sample with 
detected molecular gas. Alternatively, t he molecular gas 
withi n the galaxies could be 4) heated (jNesvadba et al.l 
IMoli . 5) kinematically prevented from collapsing int o 
GMCs (“morphological quenching”: iMartig et al.ir2009ll . 
or 6) dispersed. Here, we comment on the implications 
of our results for these latter scenarios, and on what data 
are needed to complete the picture. 

The offset observed in both Figures [6] and [7] could be 
caused by either a breakdown in the relation between CO 
and H 2 (a different value of aco): or between H 2 and the 
SFR (a different star formation efficiency), such that ei¬ 
ther is different than for the comparison galaxies. The 
burst-ending mechanisms of molecular gas heating (sce¬ 
nario 4) or morphological quenching (scenario 5) might 
alter the state of the gas relative to normal star-forming 
galaxies, acting to lower aco (see Eqn. 5) or to reduce 
the SFE for a given M{H 2 ). In the case where aco*s 
low, dispersal of the gas (scenario 6) would drive the 
post-starbursts down along the KS relation. 

Analyzing the gas state is outside the scope of this pa¬ 
per, as we cannot constrain the temperature or density of 
the CO em itting region using only the CO (1-0) and (2-1) 
lines fe.g.. iCarilli fc Walter 12013) . We need higher J„p 
lines and denser gas tracers to do so. Determining if the 



10"' 10° 10' 10^ 10^ 10'' 

^H2 0 pc ] 



6.35 5.15 3.94 2.74 1.54 0.34 

Figure 10. Top: Same as Figure [7^- (see description in ^3.3|l . 
but assuming oco = 4 for the K98 star-forming galaxies, and 
QfQO =0.8 for the K98 star bursts and our post-starburst sample. 
Here, the offset previously observed can be reconciled by assuming 
a ULIRG-type value of cx-co — 0.8 for the post-starburst sam¬ 
ple. Bottom: Same as Figure [7| but now assuming the variable 
CKCO foi* th e K98 sample and post- starburst sample using the for¬ 
mula from [Narayanan et al.l As before, triangles, circles, 

and squares represent the post-starburst sample, K98 star-forming 
galaxies, and K98 starbursts, respectively. We still observe a sys¬ 
tematic shift low on the Kennicutt-Schmidt relation, which sug¬ 
gests that the variable aco model does not resolve the differences 
between the post-starbursts and other galaxies. 

galaxy is undergoing morphological quenching requires 
resolved kinematics of the molecular gas. If the offset 
from the KS relation is not due to incorrect assumptions 
about aco (see 114.311 . the SFR 1 112.51) . or the relative spa¬ 
tial distribution of H 2 , optical light, and star formation 
f 114.21 assuming CO traces H 2 ), the intriguing possibil¬ 
ity of a lower star formation efficiency in post-starbursts 
remains. 

Another potential explanation for the offset is that 
our galaxies are still forming stars, but with a bottom- 
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heavy IMF dominated by low mass stars that are not 
detectable by our SFR indicators. Regardless of the 
IMF during the burst, any subsequent star formation 
in our sample might be expected to track that of early 
types, for which bottom-heavy IMFs have been suggested 
le.g-.. Ivan Dokkum fc Conr^ 1201(1( 1 . We test this hy¬ 
pothesis by estimating the change in the Ha-derived 
SFRs with IMF slope. For consistency, so far we have 
used the same conversion factor from L{Ha) to SFR as 
K98, wh o employed a Salpeter IMF (x = 2.35 slope). 
However, Ivan Dokkum & ConrovI (|20inl l favor a steeper, 
more bott om heavy slope of a: = 3. We use Starburst 
99 v7.0.0 (jLeitherer et al.l[20I4f) models with a variety 
of IMF slopes over the range 0.1 — IMq . For a bot¬ 
tom heavy IMF slope of a; = 3, we would underestimate 
the SFR significantly (by ~ 20x), more than enough to 
explain the offset of our galaxies from the KS relation. 

5. CONCLUSIONS 

We study molecular gas in a sample of 32 nearby 
(0.01 < z < 0.12) post-starburst (aka “E-l-A”) galaxies, 
whose optical spectra indicate a recent starburst that 
ended within the last ^Gyr. We target the CO lines 
(1-0) and (2~1) with the IRAM 30m and SMT 10m tele¬ 
scopes, constraining the molecular gas mass remaining 
in these galaxies after the starburst. Our conclusions are 
as follows: 

• Molecular gas is detected in 17 (53%) galaxies 

with CO (1-0) observations from the IRAM 30m. 
We obtain molecular gas masses of M{H 2 ) = 
1O® ®-1O®-®M0 (aco/4MQ (Kkms“^ pc^)“^), and 
molecular gas to stellar mass fractions of 
/gas ~ 10"^-10"° '^ (aco/4MQ (Kkms“^ pc^)“i), 
roughly comparable to those of star-forming galax¬ 
ies and generally larger than for early types, for a 
range of likely CO-to-H2 conversion factors (aco)- 
The upper limits on M{H 2 ) for the 15 non-detected 
galaxies range from IO^'^Mq to 10^ with the 

median more consistent with early-type galaxies 
than with star-forming galaxies. 

• We compare M{H 2 ) to the star formation rate 
(SFR), using Ha and Il„(4000) to calculate up¬ 
per limits on the current SFR in this sample. 
When compared to other star-forming, starburst¬ 
ing, and early type galaxies, the post-starbursts 
have ^ 10 — 20 X lower SFRs for a given M(H2). 

• The post-starburst sample falls ^ 4x below 
other l ocal galaxies on the Kennicutt-Schmidt re¬ 
lation (jKennicuttlll998[ l of SFR surface density vs. 
M(H2) surface density. The median locus of the 
post-starburst galaxies is offset from the relation 
defined by normal star-forming galaxies (K98) at 
5(7 significance. After considering sample selection 
effects, aperture bias, varying spatial extents of 
current star formation, optical light and H2, CO 
not tracing H2 (a different aco)i and the effect of 
IMF assumptions, we conclude the observed offset 
is likely due to suppressed SFE, a low value of aco 
consistent with ULIRGs, and/or a bottom-heavy 
IMF. 


Our results show that the end of the starburst in these 
galaxies cannot be attributed to the complete consump¬ 
tion, expulsion, or starvation of the molecular gas reser¬ 
voirs. Resolved interferometric CO maps of these galax¬ 
ies, higher J„p lines of CO for density-temperature con¬ 
straints, denser gas tracers such as HCN, and resolved 
star formation maps are necessary to more thoroughly 
study the current state of gas in post-starbursts, and to 
more accurately compare to the residual star formation. 
Understanding this possibly common phase in galaxy 
evolution will help reveal the physics of star formation in 
galaxies as well as their evolution through mergers. 
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Table 1 

Post-Starburst Targets 


Target 

R.A. 

(deg) 

decl. 

(deg) 

z 

PlatelD^ 

MJD^ 

FiberlD^ 

log(M./Mo)^ 

E>n(4000)^ 

R90I 

(arcsec) 

SFRum{Ua)2. 
M© yr-l 

SFRlirr,{Dr,(A000)) 
Mq yr-1 

ap. corr. ^ 

EAHOl 

128.640457 

17.346207 

0.0478 

2276 

53712 

444 

10.45 

1.4122 

9.35 

0.0586 

0.0149 

6.5416 

EAH02 

141.580383 

18.678055 

0.0541 

2360 

53728 

167 

9.96 

1.3798 

5.26 

0.0256 

0.0271 

3.5339 

EAH03 

222.066864 

17.551651 

0.0449 

2777 

54554 

258 

10.34 

1.4227 

7.06 

0.0196 

0.0055 

6.3128 

EAH04 

318.502258 

0.535107 

0.0269 

986 

52443 

468 

10.18 

1.2445 

8.98 

0.0674 

0.0621 

4.0061 

EAH05 

184.260117 

39.077038 

0.0653 

2001 

53493 

473 

10.00 

1.4148 

3.44 

0.0625 

0.0150 

3.6750 

EAH06 

116.456268 

31.378378 

0.0441 

755 

52235 

42 

10.53 

1.4477 

5.52 

0.3737 

0.0221 

2.8889 

EAH07 

167.824844 

11.554388 

0.0380 

1604 

53078 

161 

10.65 

1.3226 

12.50 

0.2206 

0.0292 

9.8156 

EAH08 

147.077820 

2.501155 

0.0604 

480 

51989 

580 

10.41 

1.3765 

4.48 

0.0407 

0.0176 

4.4127 

EAH09 

227.229538 

37.558273 

0.0291 

1352 

52819 

610 

10.21 

1.3668 

5.34 

0.0572 

0.0086 

4.8709 

EAHIO 

158.427979 

21.127987 

0.1053 

2376 

53770 

454 

10.24 

1.4019 

3.44 

0.0419 

0.0985 

0.9470 

EAHll 

166.419617 

5.998405 

0.0542 

1003 

52641 

87 

10.61 

1.4173 

4.67 

0.1690 

0.0271 

1.8951 

EAH12 

223.772690 

13.281012 

0.0826 

2750 

54242 

18 

10.55 

1.4002 

2.83 

0.1834 

0.0217 

1.5728 

EAH13 

155.503281 

22.163177 

0.1129 

2365 

53739 

624 

11.00 

1.2897 

4.61 

0.6212 

0.5817 

2.4663 

EAH14 

178.276855 

64.299026 

0.0622 

598 

52316 

170 

10.04 

1.2878 

3.15 

0.1005 

0.0272 

2.1201 

EAH15 

163.085205 

5.828218 

0.0411 

1001 

52670 

48 

10.40 

1.4178 

4.72 

0.0596 

0.0258 

5.6546 

EAH16 

141.740372 

42.526840 

0.1113 

870 

52325 

208 

10.74 

1.1959 

4.52 

0.5076 

5.2118 

1.5974 

EAH17 

191.215393 

-1.759901 

0.0481 

336 

51999 

469 

10.05 

1.4225 

4.67 

0.0479 

0.0634 

1.6746 

EASOl 

11.246839 

-8.889684 

0.0196 

656 

52148 

404 

10.24 

1.4919 

13.08 

0.0130 

0.0255 

5.0585 

EAS02 

49.228809 

-0.041979 

0.0231 

413 

51929 

238 

10.08 

1.4068 

11.83 

0.0253 

0.0318 

4.1478 

EAS03 

117.809624 

34.418201 

0.0628 

756 

52577 

424 

10.86 

1.4849 

7.46 

0.1710 

0.0083 

6.8937 

EAS04 

126.755821 

21.706779 

0.0153 

1927 

53321 

584 

9.99 

1.4244 

8.84 

0.0140 

0.0091 

45.2321 

EAS05 

146.112335 

4.499120 

0.0467 

570 

52266 

537 

10.57 

1.3109 

8.42 

0.1078 

0.0328 

5.4079 

EAS06 

159.488983 

46.244514 

0.0227 

962 

52620 

212 

10.14 

1.4450 

7.10 

0.0575 

0.0055 

1.3282 

EAS07 

169.781738 

58.053974 

0.0325 

951 

52398 

128 

10.54 

1.4897 

7.23 

0.0447 

0.1355 

2.7343 

EAS08 

189.900208 

12.438888 

0.0408 

1616 

53169 

71 

10.67 

1.3406 

13.55 

0.0550 

0.1720 

6.2527 

EAS09 

191.611816 

50.792061 

0.0270 

1279 

52736 

362 

10.56 

1.5223 

10.68 

0.0432 

0.0127 

4.0306 

EASIO 

196.357605 

53.591759 

0.0381 

1039 

52707 

42 

10.53 

1.5089 

6.44 

0.0189 

0.0397 

2.5848 

EASll 

242.585358 

41.854881 

0.0395 

1170 

52756 

189 

10.74 

1.5299 

10.55 

0.0765 

0.0583 

6.6460 

EAS12 

243.375778 

51.059879 

0.0336 

623 

52051 

209 

10.01 

1.1354 

18.20 

0.0314 

0.2572 

1.3772 

EAS13 

246.760666 

43.476093 

0.0462 

815 

52374 

586 

10.95 

1.4548 

11.41 

0.0860 

0.0101 

7.2333 

EAS14 

316.286133 

-5.399832 

0.0826 

637 

52174 

584 

11.31 

1.5036 

9.68 

0.4120 

0.0396 

4.3280 

EAS15 

343.778320 

0.977756 

0.0533 

379 

51789 

579 

10.83 

1.3498 

9.70 

0.1697 

0.5866 

4.0352 


JNotes: 


^SDSS Spectra Identification from DR7 dAbazaiian et ani2009l‘) . 

'^MPA-JHU Stellar Masses, method described in IBrinchmann et al.l H2004l 'l. 

1 ) 77 .( 4000 ) from SDSS spectra. 

^Petrosian 90% size in r band, from SDSS photometry. 

^Limit on SFR from Ha measurements, before aperture correction. Conv ersion from 
IKennicutt et aO l| 19941') . line fluxes from MPA-JHU catalogs l|Aihara et aLlf^Ollh . See ^2.41 
^Lim it on SFR from Z)n(4000) measu reme nts, from MPA-JHU catalogs HBrinchmann et al.l 
(2001) . before aperture correction. See 3221 

^Aperture correction to conve rt fib er-based SFRs to global SFRs, from MPA-JHU catalogs 
l|Brinchmann et al.ll200^ . See ^2.51 
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Table 2 

IRAM 30m CO (1-0) observations 


Target ^ 

^obs 

(hours) 

Icol 

(K km s“^) 

(10^ Kkms“^ pc^) 

M{R2]^ 

(10" Mb ) 

FWHM_^ 

(km s~^) 

EAHOl 

0.30 

2.58± 0.33 

128.6± 16.4 

514.5± 65.7 

309.8± 18.4 

EAH02 

0.31 

1.32± 0.27 

84.3± 17.4 

337.0± 69.7 

172.9± 29.8 

EAH03 

0.20 

3.60± 0.39 

158.0± 17.2 

632.1± 68.9 

158.7± 8.6 

EAH04 

0.51 

0.59± 0.13 

9.2± 2.0 

36.7± 7.9 

106.7± 20.3 

EAH05 

0.60 

0.97± 0.21 

90.7± 19.5 

362.7± 77.9 

278.7± 44.8 

EAH06 

0.71 

< 0.59 

< 24.84 

< 99.34 


EAH07 

1.02 

< 0.33 

< 10.31 

< 41.25 


EAH08 

0.83 

0.45± 0.15 

35.6± 12.0 

142.5± 47.8 

233.7± 56.9 

EAH09 

0.61 

0.43± 0.13 

7.8± 2.4 

31.4± 9.7 

119.5± 28.5 

EAHIO 

0.43 

0.73± 0.20 

181.5± 48.6 

726.1± 194.3 

275.1± 58.8 

EAHll 

0.40 

< 0.60 

< 38.53 

< 154.12 


EAH12 

0.70 

< 0.44 

< 65.87 

< 263.47 


EAH13 

0.94 

0.68± 0.13 

193.8± 37.6 

775.2± 150.3 

190.1± 27.8 

EAH14 

0.80 

< 0.53 

< 44.96 

< 179.85 


EAH15 

0.60 

< 0.75 

< 27.57 

< 110.27 


EAH16 

0.30 

< 0.56 

< 155.85 

< 623.40 


EAH17 

0.43 

< 0.51 

< 25.67 

< 102.66 


EASOl 

0.33 

< 0.76 

< 6.23 

< 24.93 


EAS02 

0.40 

l.llit 0.34 

12.8± 3.9 

51.2± 15.5 

162.7± 44.0 

EAS03 

0.47 

1.66± 0.24 

143.6± 21.1 

574.4± 84.5 

270.5± 30.3 

EAS04 

1.56 

< 0.28 

< 1.38 

< 5.51 


EAS05 

0.74 

0.64± 0.20 

30.5± 9.5 

121.8± 38.2 

346.9± 82.0 

EAS06 

0.31 

3.83± 0.39 

42.6± 4.3 

170.6± 17.2 

112.9± 4.0 

EAS07 

0.61 

< 0.47 

< 10.82 

< 43.27 


EAS08 

1.14 

< 0.28 

< 10.06 

< 40.24 


EAS09 

0.54 

2.12± 0.27 

33.3± 4.3 

133.3± 17.2 

267.8± 20.5 

EASIO 

0.80 

< 0.49 

< 15.53 

< 62.12 


EASll 

0.40 

< 0.51 

< 17.20 

< 68.80 


EAS12 

1.10 

0.36± 0.12 

8.7± 2.8 

34.9± 11.3 

141.7± 29.1 

EAS13 

0.39 

< 0.70 

< 32.75 

< 131.01 


EAS14 

1.61 

0.83± 0.17 

124.6± 25.3 

498.5± 101.1 

428.7± 73.4 

EAS15 

0.74 

0.48± 0.15 

29.9± 9.5 

119.8± 38.0 

166.2± 57.2 


Notes: 

^Lines in bold represent > 3 (j detections in IRAM 30m CO (1-0) observations. 
^Upper limits are shown at the 3cr level. 

'^Calculated using Lqq = 23.5 12^ Ico (1 + z)~^. 

'^Masses calculated assuming (Kkms“^ pc^)“^, M(H 2 ) = 

®FWHM from Gaussian fit to data. 
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Table 3 

IRAM 30m and SMT CO (2-1) observations 


Target 

tobs^ 

(hours) 

-L b 

^IRAM _ 

(K km s“^) 

IRAM 30m FWHM_^ 
(km s“^) 

4 , d 
^OOS ^ 

(hours) 

-1 b 

(K km s“^) 

SMT EWHM_^ 
(km s“^) 

(arcsec) 

EAHOl 

0.30 

4.67±0.57 

337.7±18.4 

3.10 

1.12±0.27 

275.7±31.2 

14.6^4? 

EAH02 

0.26 

1.39±0.46 

81.6±18.2 





EAH03 

0.20 

5.60±0.61 

144.0±6.3 

2.85 

0.64±0.15 

149.5±25.2 

SAtli 

EAH04 

0.41 

1.51±0.23 

125.3±17.0 

6.60 

0.20±0.03 

126.5±16.9 


EAH05 

0.54 

1.21±0.30 

266.4±54.0 

6.00 

<0.72 



EAH06 

0.71 

<1.53 






EAH07 

1.03 

<0.48 






EAH08 

0.76 

<1.06 


3.80 

<0.18 



EAH09 

0.61 

0.78±0.18 

127.0±24.2 

9.50 

<0.20 



EAHIO 

0.43 

<1.83 






EAHll 

0.40 

<0.66 






EAH12 

0.70 

<1.15 






EAH13 

0.94 

<1.17 






EAH14 

0.80 

<0.98 






EAH15 

0.60 

<1.14 






EAH16 

0.32 

<2.05 






EAH17 

0.43 

<1.69 






EASOl 

0.34 

<0.81 






EAS02 

0.40 

2.31±0.70 

202.2±62.5 

4.60 

0.38±0.09 

131.2±16.9 


EAS03 

0.47 

1.84±0.43 

71.0±18.0 

6.80 

<0.28 



EAS04 

1.61 

<0.47 






EAS05 

0.75 

<0.81 


4.60 

<0.77 



EAS06 

0.31 

5.19±0.54 

117.6±3.6 

1.50 

2.10±0.21 

111.7±8.0 


EAS07 

0.61 

<0.74 






EAS08 

1.00 

<0.61 






EAS09 

0.55 

3.39±0.39 

276.1±13.8 

2.60 

0.79±0.26 

234.6±46.3 


EASIO 

0.80 

<0.82 






EASll 

0.39 

<0.52 






EAS12 

1.10 

0.62±0.14 


5.55 

<0.27 



EAS13 

0.44 

<0.81 






EAS14 

1.61 

1.15±0.30 

598.1±199.7 





EAS15 

0.70 

1.04±0.22 

299.0±58.1 

4.50 

<0.16 




Notes: 

^Time on source (hours) IRAM 30m. 

^Upper limits are shown at the 3cr level. 

‘^FWHM from Gaussian fit to data. 

‘^Time on source (hours) SMT. 

^Approximate Gaussian source size (FWHM) of GO emitting region, see text f H4.2ll for 
details 
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Figure 11. CO(l-O) and CO(2-l) spectra from IRAM and SMT for galaxies with IRAM-30m CO (1-0) detections (> 3cr) in our post- 
starburst sample. Spectra are shown in units of both main beam temperature [Jy]- Grey lines show the unbinned IRAM 

data for 5 km/s channels, and black lines show the data binned to 20 km/s. Dashed red lines represent the rms of the binned data. SMT 
data are shown in 13 km/s bins. Blue horizontal lines at bottom represent the integration intervals, as described in the text. 
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Figure 11. continued 
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Figure 11. continued 
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Figure 11. continued 
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Figure 12. CO(l-O) and CO(2-l) spectra from IRAM and SMT for galaxies with CO (1-0) not detected with the IRAM 30m. Spectra 
are shown in units of both main beam temperature [rifiK] and Si/ [Jy]. Grey lines show the unbinned IRAM data for 5 km/s channels, 
and black lines show the data binned to 20 km/s. Dashed red lines represent the rms of the binned data. SMT data are shown in 13 km/s 
bins. Blue horizontal lines at bottom represent the integration intervals, as described in the text. 
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Figure 12. continued 
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Figure 12. continued 











































































































